Observable NMR signal from circulating current order in YBCO 
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Assuming, as suggested by recent neutron scattering experiments, that a broken symmetry state with orbital 
current order occurs in the pseudo-gap phase of the cuprate superconductors, we show that there must be asso- 
ciated equilibrium magnetic fields at various atomic sites in the unit cell, which should be detectable by NMR 
experiments. 



I. INTRODUCTION 

Varma 1 has proposed that the pseudogap state of the high 
temperature superconducting cuprates is a broken symmetry 
phase with equilibrium circulating currents within the unit cell 
as shown schematically in Figure [T] Indeed, recent polarized 
neutron scattering data 2 3 suggest that in at least some of the 
cuprates there is a phase transition at a pseudo-gap tempera- 
ture, T* > T c , to a state with a form of magnetic order which 
preserves the translational symmetries of the crystal. (Here 
T c is the superconducting transition temperature and T*, de- 
termined from neutron scattering, is comparable to the previ- 
ously determined pseudo-gap crossover temperature derived 
from transport and NMR studies.) 

Specifically, the neutron scattering data suggest that there 
is some form of intra-unit-cell antiferromagnetic ordering in 
the pseudo-gap phase of underdoped YBCO and Hgl201 with 
ordered moments with magnitudes of order . 1 /is , tilted away 
from the c axis by a substantial angle (roughly 35° — 65°). 
This is broadly compatible with the Varma loop order or with 
an alternative model (of a sort proposed by Fauque et a/pJwith 
ordered moments on O sites. 

Here, we determine the implications for Cu, Ba, and O 
NMR of the existence of magnetic symmetry breaking of the 
sort suggested by the neutron scattering experiments. Specif- 
ically, if we assume Varma loop order and take the sim- 
plest physically plausible model of the distribution of currents 
within the unit cell, we can analyze the neutron data to obtain 
quantitative estimates of the magnitudes of the putative orbital 
currents. From these, we can compute the expected value of 
the associated magnetic fields at various sites in the unit cell. 
The resulting fields are large enough to be detectable, so this 
constitutes a falsifiable consequence of the model. 



II. SYMMETRY CONSIDERATIONS 

We begin by briefly reviewing the symmetry considerations 
governing the existence of a static (thermodynamic) magnetic 
field, B(r) at any given position in the unit cell of a crystal. 
Specifically, if certain symmetries are preserved, it is possible 
to prove that the magnetic field vanishes. If the field does not 
vanish by symmetry, the likely implication is that the field is 
non-zero. 

Let Qf be the set of point group symmetry operations in the 
considered state of the system which leave the spatial point f 



invariant. If there exists any group element g £ Q?, such that 
g transforms B a (r) — > —B a (r), it follows that B a (r) = 0. 

Since B is odd under time-reversal (T), if time-reversal 
symmetry is unbroken, then B(f) = for all r. However, 
if time-reversal symmetry is broken, the field may still van- 
ish at points of sufficiently high symmetry. The Varma state 
breaks both time-reversal and rotation by tt (R^) about a line 
parallel to the c axis and passing through any atomic site, but 
preserves the product of these (TR^). Since TR^ transforms 
B c — > — B c , the out-of-plane component of the magnetic field 
must vanish at all atomic sites for both Hg-1201 and YBCO. 

Because B is a pseudo-vector, any component of B that lies 
in a mirror plane is odd under reflection through this plane. 
Thus, if there is a symmetry with respect to a reflection plane 
passing through point r, then all components of B(f) within 
that plane must vanish. We will assume that the loop order 
leaves unbroken the reflection symmetry about the Cu-O plane 
in Hg-1201, and about the Y plane in YBCO, both of which 
we term M. Along with TR„ symmetry, this assures that 
there is no net ferromagnetism in the loop ordered state. This 
reflection symmetry also implies that the in-plane components 
of the magnetic field vanish in the Cu-O plane of Hg-1201. 

Finally, to simplify the discussion, we will assume the exis- 
tence of an additional unbroken reflection symmetry D about 
a diagonal plane containing the Cu, apical O, and Ba sites in 
both materials. This is at best approximate in YBCO, which 
is orthorhombic in the normal state, and thus would be mon- 
oclinic in the Varma state. However, since the YBCO lattice 
is tetragonal to within 2%, deviations from this symmetry are 
likely to be small. 

In summary, if there exists a Varma loop ordered state in 
Hg-1201, the M and TR^ symmetries require B{r) to vanish 
at the Cu and planar O sites. However, non-vanishing fields 
parallel to the Cu-O plane are permitted by symmetry at the 
apical O and Ba sites of Hgl201, and at the Ba, Cu, apical 
O, and planar O sites of YBCO. The D reflection symmetry 
further constrains the directions of the fields. For instance, 
in Hg-1201, the field at the apical O and Ba sites must point 
along a unit cell diagonal. (Note that similar symmetry con- 
siderations applied to the d-density wave state 5 imply vanish- 
ing magnetic fields at the Cu and all the O sites in Hg-1201 
and at the Cu, and apical O sites in YBCO, but permit an in- 
plane magnetic field at the in-plane O sites in YBCO, and an 
out-of-plane magnetic field at the Ba sites of both Hg-1201 
and YBCO.) To estimate of the order of magnitude of the ex- 
pected fields, we consider the field at the center of a circular 
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current loop of dipole moment /i and radius 7-: 

where a — 3.8Ais the Cu-Cu distance. Taking r = | and /_i = 
0.1/zb (as reported by Mook et al®), this gives B ~ 250G. 

HI. EXPLICIT MODEL 

For simplicity, we consider a model in which only the Cu- 
O planes are electronically active, and in which we ignore 
the orthorhombicity of the YBCO lattice and any small buck- 
ling of the in-plane O-Cu-O bond. We assume the symmetries 
mentioned in the previous section as well as lattice translation 
symmetry. 

The ideal two dimensional cartoon of the state originally 
proposed by Varma is sketched in Figure [T] Since in the 
cuprates, the distance to the apical O is comparable to the dis- 
tance between planar O's, we must take into account currents 
involving the apical O, even if we take the simplest (shortest- 
range) version of this state adapted to the actual materials. 
(This was recognized previously in RefspEla as it is neces- 
sary to account for the presence of the in-plane components 
of the magnetic order inferred from the neutron scattering ex- 
periments.) The assumption of near-neighbor currents and the 
symmetries previously assumed lead to the pattern of currents 
shown in Figures [2] and [3] for YBCO and Hg 1201, respec- 
tively. Note that TR n symmetry forbids a current between 
the apical O and Cu sites, while the reflection symmetries M 
and D ensure that the precise pattern of broken symmetry in 
the proposed state is defined by the three independent currents 
labeled Ij in the figures. 

The model presented is broadly consistent with the neutron 
scattering data in that it would result in magnetic scattering 
intensity at suitable Bragg vectors with a polarization depen- 
dence reflecting comparable strengths of the in-plane and out- 
of -plane magnetic fields. However, there is at least one un- 
resolved discrepancy: The in-plane component of the mag- 
netization at (100) must vanish so long as M is unbroken. 
Experimental data 2 find a tilted moment at (100), though this 
is only reported for one experiment in the published literature, 
and with a substantial margin of error. 

Before making a quantitative comparison with experiment, 
we have imposed the additional constraint, not required by 
symmetry, that there be no net current flowing through the sys- 
tem in equilibrium. This is equivalent to taking I3 = I\ + I2 
in Figures [2] and [3] Assuming that "form factor" effects as- 
sociated with the spatial extent of these currents can be ne- 
glected, the resulting magnetic field that would be produced 
by this state can be directly computed (as sketched in the Ap- 
pendix) in terms of the two remaining independent parame- 
ters Ji and I 2 - These can be determined by comparing the 
measured and predicted spin-flip magnetic scattering cross- 
sections for orthogonal neutron polarizations at a single Bragg 
wave-vector. As the strongest scattering (and consequently 
the smallest experimental uncertainty) occurs at (011) in both 




Figure 1 : Top view of circulating currents in the Varma loop ordered 
state. Cu sites are black rectangles, O sites white circles. 




Figure 2: Microscopic model for YBCO. O sites are spheres, Cu sites 
are cubes. 



YBCO and Hg-1201, we determine the values of Ij for each 
material using data from that peak only. Since the scattering 
cross-section is quadratic in the currents, there are two inde- 
pendent solutions for Ii and I 2 - Naturally, once the currents 
are determined, it is straightforward to compute the field at 
any particular spatial point. We quote results for each of the 
two independent solutions for 1\ and I 2 . 



IV. RESULTS AND DISCUSSION 

For the case studied in RefPof YB 2 C306+«5 with 5 = 0.6 
and T c = 61K, the analysis described above leads to the pre- 
diction of fields of 100's of G at various sites in the unit cell, 
consistent with the dimensional estimate presented above. In 
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Figure 3: Microscopic model for Hgl201. O sites are spheres, Cu 
sites are cubes. 



particular, we predict fields of 700G or 710G at the copper 
site, 230G or 130G at the barium site, 170G or 280G at the 
apical oxygen sites, and 360G or 350G at the in-plane oxygen 
sites. Additional symmetries in Hgl201 lead to quite differ- 
ent results. In particular, as required by symmetry, the fields 
at both the Cu and the planar O sites vanish exactly. How- 
ever, if we apply the present analysis to the neutron data of 
RefPfor underdoped Hg 1201 (T c = 81K), we infer that at 
low temperatures there should be a field of 180G or 200G at 
the apical oxygen site, and 240G or 170G at the Barium site. 
Observed^ Cu NMR line-widths in YBCO are of order 100G 
and O linewidths substantially smaller, so these effects should 
be readily observable. 

There are, of course, uncertainties in the quantitative esti- 
mates we have made associated with the error bars on the neu- 
tron data. Most importantly, the assumed pattern of currents is 
not fully consistent with the neutron data; as mentioned above, 
it follows from the assumed symmetries of the state that the 
in-plane moment at (100) must vanish, whereas a small, in- 
plane component is reported in experiment 2 . There are also 
theoretical uncertainties associated with the simplifying as- 
sumptions that we have made in defining the model used for 
explicit calculations. Specifically, some quantitative changes 
could arise from form factors associated with the non-zero 
spatial extent of the equilibrium currents, and with currents 
flowing between further neighbor sites in the lattice. We do 
not believe, but have not proven, that these uncertainties will 
not cause qualitative revisions in the quantitative estimates 
we have made. We have also treated the currents classically, 
which is a valid procedure in typical situations involving bro- 
ken symmetry states; however, Varma 8 has argued that there 
is an essentially quantum character to the loop ordered state 
which invalidates such an approach. Finally, a magnetic order- 
ing that appears static within the frequency resolution of neu- 
tron scattering may be rapidly fluctuating on NMR timescales. 
Though no obvious mechanism exists for an order fluctuating 



on intermediate timescales, it would be consistent with the 
neutron results and the symmetries assumed, yet yield a null 
result in NMR experiments. In any case, the presence of static 
magnetic fields at the various atomic sites in the crystals, if 
seen, would constitute strong evidence of the putative current 
loop order in the pseudo-gap phase of these materials. 

At present, we are unable to find published Cu or O NMR 
studies which give clear results concerning the existence or 
absence of fields of the predicted magnitude in YBCO or 
Hgl201. Work by Strassle et al. 9 sets an upper bound of less 
than 1G for the field at the Barium site of YBCO 248. How- 
ever, since neutron scattering is not available for this mate- 
rial, this bound cannot be directly compared with a predicted 
magnitude. We hope that the results in this Brief Report will 
encourage experimental tests of local magnetic fields in the 
pseudogap state of the aforementioned materials. 
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V. APPENDIX: SKETCH OF CALCULATIONS 

The neutron scattering data used in these calculations are 
the cross sections for elastic magnetic neutron scattering in the 
spin-flip channel for two independent polarizations. The Born 
approximation gives the differential scattering cross section at 
momentum transfer K (a reciprocal lattice vector) as 

= ©)'(™)>*-*(*- 3 *)f 

where g = —3.826 is the neutron g-factor, /.in the nuclear 
magneton, P the polarization direction of the incoming neu- 
trons, and Bg is the Fourier transform of the magnetic field. 

The momentum transfer K employed in Refers (0,1,1). 
The magnetic field component at this wavevector can be ex- 
pressed in terms of the currents I\, 1%, and 1% as 

^ Ani K x 



where if is a reciprocal lattice vector, and the current density 
Jg can be expressed as 

sin ( K ■ e„ ) 

n A ' e " 

In this expression the sum runs over the current-carrying 
bonds (segments) in a unit cell. For a given segment n, Ij n 
is the current (J\, I2, or ^3, as defined previously), e n is the 
direction of current flow, L n is the length of the bond, and 7Z n 
is the position of the bond midpoint. 

With these results, we can write the scattering cross section 
as a function of the currents, and then invert to solve for these 
currents in terms of the experimental data. We then evaluate 
the real space magnetic field using the Biot Savart Law. All 
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told, the field is 

B(r) 



a x e r 




In this expression the first sum runs over all unit cells of the 
crystal (i.e. R runs over all Bravais lattice vectors), the second 
sum is over the current-carrying bonds in a unit cell as above, 
and the position of a given segment is a = a n ^ = R—r+lZ n . 
The sum converges rapidly. 
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